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Comparison  of  Soil  Freezing  Curve  and  Soil  Water  Curve  Data 
for  Windsor  Sandy  Loam 


PATRICK  B.  BLACK  AND  ALLEN  R.  TICE 


INTRODUCTION 

The  CRREL  Outdoor  Land  Treatment  Research 
Facility  consists  of  six  test  beds  containing  two  soil 
types.  The  soils’  ice-free  physical  and  chemical 
properties  have  been  reported  by  Iskandar  et  al. 
(1979).  Our  report  presents  recently  obtained  lab¬ 
oratory-measured  unfrozen  water  content  data  as 
a  function  of  temperature,  plotted  as  soil  freezing 
curves  (SFC),  for  all  three  horizons  of  the  Wind¬ 
sor  sandy  loam  present  in  the  test  beds.  It  then 
demonstrates  how  such  a  soil  data  base  is  related 
to  the  previously  measured  ice-free  water  content 
data  as  a  function  of  matric  suction,  plotted  as  soil 
water  curves  (SWC).  The  theoretical  basis  for  the 
similarity  of  the  SFC  and  SWC  is  understood 
(Miller  1965,  Koopmans  and  Miller  1966)  but  the 
value  of  relating  different  data,  collected  by  dif¬ 
ferent  researchers  at  different  times  for  the  same 
field  soil,  needs  to  be  shown  to  illustrate  the  use¬ 
fulness  of  soil  data  bases. 

This  report  first  presents  a  generalized  notation 
for  the  state  of  soil  water  and  then  reviews  the  sim¬ 
ilarity  between  the  state  of  liquid  water  in  unfro¬ 
zen  soil  and  air-free  frozen  soil  in  the  generalized 
notation.  Unfrozen  water  content  data  of  the 
Windsor  sandy  loam  measured  by  nuclear  mag¬ 
netic  resonance  (NMR)  are  compared  to  data  on 
water  content  as  a  function  of  matric  pressure  pre¬ 
sented  by  Iskandar  et  al.  (1979)  through  the  simi¬ 
larity  relationships  presented. 


SOIL  VARIABLE  <f> 

For  discussion  of  the  state  of  water  in  soil-water 
systems,  it  is  convenient  to  introduce  the  variable 
4>jk: 

*jk  =  «j-"k.  0) 


the  pressure  difference  between  water  in  the  two 
phases  j  and  k  (Black  and  Miller  1985,  Black 
1985).  The  matric  pressure  for  ice-free  soil  is  ex¬ 
pressed  as 

^aw  =  ua  —  uw  (2) 

where  the  subscripts  a  and  w  refer  to  soil  air  and 
water  pressures,  respectively.  Likewise,  the  state 
of  water  in  an  air-free  frozen  soil  is  expressed  as 

$  iw  —  U[  —  Mw  (3) 

where  «j  is  the  ice  pressure.  Most  often,  though, 
the  state  of  water  in  air-free  frozen  soil  is  expressed 
in  terms  of  temperature,  0  (in  Celsius  degrees).  If 
the  soil  is  devoid  of  solutes,  the  Clapeyron  equa¬ 
tion  can  be  expressed  as: 


where  y)  is  the  specific  gravity  of  ice  and  h  is  volu¬ 
metric  latent  heat  of  fusion.  Solving  eq  3  and  4  for 
Uj  gives  the  connection  between  water  pressure  uw, 
ice  temperature  0  and  <t> jw: 

*iw  =  (7i-l)«w  (5) 

Another  useful  classification  to  employ  in  dis¬ 
cussing  the  physical  behavior  of  soil  includes  the 
concepts  of  adsorption  space  and  capillary  space. 
Adsorption  space  is  that  zone  in  which  soil  water 
is  strongly  affected  by  surface  forces  (real  or  vir¬ 
tual)  emanating  from  the  soil.  Capillary  space  is 
the  remaining  zone  in  which  soil  water  is  not  af¬ 
fected  by  soil  force  fields  but  is  governed  by  the 
laws  of  surface  tension.  Granular  soils  contain 
mostly  capillary  water,  while  highly  colloidal  soils 


are  dominated  by  adsorption  forces.  The  nature 
and  formulation  of  the  adsorption  space  are  not 
addressed  in  this  report  but  are  employed  only  as  a 
general  classification  scheme  for  soils. 

SWC  AND  SFC  SIMILARITY 

Miller  (1965)  hypothesized  that  if  the  same 
states  of  soil  moisture  content  and  distribution  are 
achieved  by  a  freezing  and  thawing  process  as  in  a 
drying  and  wetting  process,  then  the  two  states 
should  be  similar  and  interchangeable  if  the  soil  is 
either  colloidal  (adsorption  space  >>  capillary 
space)  or  colloid-free  (adsorption  space  <  <  capil¬ 
lary  space).  In  the  former  case,  the  ice-free  SWC 
data  and  air-free  SFC  data  for  the  same  soil  at  the 
same  density  and  subjected  to  similar  histories 
should  be  directly  related  (i.e.  </>aw  =  0 jw).  In  the 
latter  case,  an  additional  correction  factor  that 
takes  account  of  the  differences  in  surface  tension 
(oraw,  <7jw)  is  needed  to  bring  the  data  into  coinci¬ 
dence  [i.e.  4>aw  =  (uaw/ O|W)0jW]. 

The  experiments  of  Koopmans  and  Miller 
(1966)  prove  this  hypothesis  when  the  same  soil  is 
used  to  measure  SWC  and  SFC  data.  They  found 
that  the  transformation  from  temperature  to  <ftw 
given  by  eq  5  was  sufficient  to  directly  relate  SFC 
and  SWC  data  for  soils  dominated  by  adsorption 
forces.  In  soils  dominated  by  capillary  space,  their 
empirically  determined  ratio  of  air-water  to  ice- 
water  surface  tensions  was 

^  =  -2=-  =  2.20.  (6) 
•Piw  ffiw 

Table  1.  Relationships  between  <£aw  and  01.  for 
similar  Hqnid  water  contents. 

Adsorption  space  «  capillary  space 

0IW  ~  (Tl) 

or 

kPa 

<t>„  (kPa)  =  (2.2)  -1110  —  6  (°C)  <T2) 

Adsorption  space  »  capillary  space 

6*m  -  •hm  f^) 

or 

0„(kP*)  =  -1110^  6  ("Q  (T4) 


With  a  standard  textbook  value  for  <raw  of  72.7  x 
10'*  J/m\  their  predicted  </jw  of  33.1  x  10'*  J/m1 
agreed  fairly  well  with  estimates  obtained  indepen¬ 
dently  of  soil  physics  investigations  (Hesstvedt 
1964). 

Equations  5  and  6  allow  the  interchange  of  ice- 
free  SWC  and  air-free  SFC  data  for  soils  that  are 
dominated  by  either  absorption  space  or  by  capil¬ 
lary  space.  Table  1  contains  the  relationships  be¬ 
tween  <>aw  and  4>\*,  by  assuming  uw  to  be  zero 
(gauge)  and  choosing  values  for  the  physical  cons¬ 
tants  in  eq  5  for  the  two  soil  classifications. 

MATHEMATICAL  REPRESENTATION  OF 
SWC  AND  SFC  DATA 

Brooks  and  Corey  (1964)  proposed  a  function 
to  represent  SWC  data;  from  a  large  number  of 
observations  for  ice-free  soils,  they  found  that  the 
relative  degree  of  saturation  or  dimensionless 
water  content 

IF-  Wd 

s  =  ivfd  <7> 

could  be  reasonably  described  by  the  relationship: 

s -(•£)“  <8> 

where  IF,  Wd  and  IF,  are,  respectively,  the  water 
contents  (usually  volumetric)  at  a  given  4>aw,  at  the 
lower  limit  of  drying  and  at  saturation.  <t>\,  is  the 
air  entry  value  and  6  is  a  free  parameter  determined 
from  a  "curve  fit”  to  the  data  and  appears  to  be 
related  to  the  pore-size  distribution  of  the  soil. 
Small  values  of  S  are  found  to  correspond  to  soils 
with  a  wide  range  of  pore  sizes;  large  6  values  are 
obtained  when  grain  sizes  are  nearly  uniform. 
Equation  8  is  valid  for  the  range,  </>tw  >  fa  other¬ 
wise,  S  =  1  (i.e.  W  -  for  4>aw  <  fa 

The  predictive  accuracy  of  eq  8  depends  upon 
the  value  for  1Yd,  the  lower  limit  of  drying.  Unfor¬ 
tunately,  data  at  very  large  values  of  4>aw  are  sel¬ 
dom  collected  because  of  experimental  complica¬ 
tions  requiring  a  value  for  Wd  to  be  inferred  by  ex¬ 
tending  outside  the  range  of  data  through  numeri¬ 
cal  or  graphical  procedures.  When  optimal  values 
for  each  parameter  are  determined  by  nonlinear 
optimization  techniques  for  experimental  data, 
negative  values  for  IFd  are  obtained  (Fields  et  al. 
1984).  Graphical  interpolation  methods  (Brooks 
and  Corey  1964,  van  Genuchten  1978)  or  rule  of 
thumb  guesses  are  biased  at  best. 
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Equation  8  can  also  be  employed  to  represent 
SFC  data  (Black  1985,  Black  and  Miller  1985)  but 
the  power  curve  relationship  proposed  by  Ander¬ 
son  and  Tice  (1973)  is  most  often  used.  Their  rela¬ 
tionship  between  gravimetric  water  content  and 
the  absolute  value  of  temperature  ( |  °C  | )  can  be 
written  more  generally  as  (on  a  volumetric  or  grav¬ 
imetric  basis) 


W=  a 


(9) 


where  a  and  0  are  free  parameters  determined 
from  a  log-log  regression  of  the  data  and  the  equa¬ 
tions  in  Table  1  are  used  to  relate  temperature  and 

^ iw 

Implicit  in  eq  9  is  the  assumption  that  the  mini¬ 
mum  value  of  water  content  is  zero  (i.e.  Wd  =  0) 
at  an  infinitely  cold  temperature.  This  would  seem 
to  be  reasonable  and  is  corroborated  by  the  suc¬ 
cessful  fit  of  data  to  the  equation  (see  App.  A). 
Similarly,  it  would  seem  reasonable  that  Wd 
should  also  be  zero  at  infinite  <£aw,  but  the  relative¬ 
ly  small  magnitude  (several  bars)  encountered  in 
determining  SWC  data  leaves  the  impression  that 
there  is  finite  water  content  even  at  very  large  0aw. 
When  SFC  data  are  included  with  SWC  data,  as  is 
done  below,  eq  9  is  sufficient  to  describe  all  the 
data  greater  than  the  air-entry  value. 

When  WA  is  assumed  to  be  zero,  eq  8  and  9  are 
equivalent  and  the  free  parameters  in  eq  9  take  on 
physical  significance.  In  this  case: 


0  =  -5;  a 


;  when  JFd  =  0 


The  air-entry  value  can  be  obtained  from  the  sim¬ 
ple  log-log  regression  of  data  if  W,  is  known. 


tect  only  hydrogen,  was  used  to  determine  unfro¬ 
zen  water  contents.  Each  test  specimen  was  re¬ 
moved  sequentially  from  the  bath,  wiped  dry  and 
inserted  in  the  NMR  probe.  After  the  four  seconds 
required  to  record  specimen  temperature  and  first 
pulse  NMR  signal  amplitude,  the  specimen  was  re¬ 
inserted  in  the  bath.  When  all  observations  had 
been  completed  at  a  given  temperature,  the  bath 
temperature  was  decreased.  Thermal  equilibrium 
was  attained  after  approximately  45  minutes  and 
measurements  repeated. 

When  the  bath  temperature  reached  -25°C,  the 
cooling  was  stopped  and  a  warming  cycle  started. 
Warming  observations  were  made  in  a  manner 
similar  to  those  determined  during  cooling,  and 
measurements  stopped  when  specimens  were  com¬ 
pletely  thawed.  Water  contents  were  then  deter¬ 
mined  gravimetrically  and  are  presented  in  Table  2 
along  with  dry  bulk  densities,  volumetric  water 
contents,  porosities  and  relative  degrees  of  satura¬ 
tion  based  upon  the  8-cm’  specimen  volume. 

The  unfrozen  water  contents  at  each  tempera¬ 
ture  were  determined  from  the  measured  first 
pulse  NMR  amplitude.  The  technique  is  described 
by  Tice  et  al.  (1978,  1981,  1982).  Briefly,  the  ratio 
of  the  gravimetric  water  content  to  the  first  pulse 
amplitude  of  the  ice-free  case  for  each  test 
specimen  was  determined.  Unfrozen  water  con¬ 
tents  were  then  deduced  by  multiplying  the  meas¬ 
ured  first  pulse  amplitude  at  the  different  temper¬ 
atures  by  the  above-determined  ice-free  ratio.  The 
data  and  regression  analysis  for  eq  9  are  listed  in 
Appendix  A. 

CHARACTERIZATION  OF  SWC 

Volumetric  water  content  ( fVy)  at  various  matric 
suctions,  4>aw,  for  the  three  horizons  of  the  Wind- 


NMR  MEASUREMENT  OF 
UNFROZEN  WATER  CONTENT 

Soil  samples  from  all  three  horizons  of  the 
Windsor  soil  were  collected.  Remolded  specimens 
were  prepared  to  different  gravimetric  water  con¬ 
tents  and  packed  to  different  bulk  densities  in 
sealed  plastic  tubes  at  a  volume  of  8  cm1.  These 
specimens  were  then  placed  in  a  constant  tempera¬ 
ture  bath  at  approximately  ~6°C  to  induce  freez¬ 
ing,  at  which  time  the  bath  temperature  was  in¬ 
creased  to  a  temperature  just  below  0°C  and  al¬ 
lowed  to  equilibrate  overnight. 

A  Praxis  model  PR-103  pulsed  nuclear  magnetic 
resonance  (NMR)  analyzer,  factory-tuned  to  de¬ 


TaMe  2.  Initial  physical  characteristics  for  Wind¬ 
sor  soil  determined  by  Iskander  et  al.  (1979). 

Soil 


horizon 

w»r 

8</ 

w,r 

o. 

c 

S 

A 

32.7 

1.376 

45.0 

2.63 

47.7 

94.3 

B 

22.1 

1.582 

35.0 

2.69 

41.2 

85.0 

C 

23.3 

1.334 

36.0 

2.73 

43.8 

82.2 

WwT  -  total  gravimetric  water  content  (mass  water /mass  soil) 
x  100 

m  total  volumetric  water  content  (volume  water/total 
volume)  x  100 

ed  *  dry  density  (mass  soil/total  volume) 

C,  m  specific  gravity 

e  -  porosity  (volume  pores/total  volume)  x  100 
S  «  relative  degree  of  saturation  (K,/e)  x  100 
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sor  soil  were  obtained  from  Figures  2,  3  and  4  of 
Iskandar  et  ai.  (1979)  and  listed  in  Appendix  C. 
The  initial  physical  characteristics  are  listed  in 
Table  2. 

DISCUSSION 

Ice-free  SWC  and  air-free  SFC  data  can  be  com¬ 
pared  only  on  the  basis  of  the  same  soil  at  the 
same  bulk  density  and  with  similar  histories  for 
each  (i.e.  freezing  and  drying  curves  or  wetting 
and  wanning  curves  can  be  compared).  In  more 
general  terms,  the  soil  must  meet  similar  media  re¬ 
quirements  (Miller  and  Miller  19S6).  Since  the 
same  soil  was  used  for  both  SFC  and  SWC  meas¬ 
urements,  only  specimens  with  identical  bulk  den- 


Tabie  3.  Summary  statistics  of  combined  SFC 
and  SWC  data  in  terms  of  If,  (W)  and 
(kPa)  fitted  to  eq  9  (i.e.  W,  = 


Soil 

horizon 

a 

0 

r' 

*6 

OtPa) 

A 

D/C 

W/w 

14.709 

12.021 

-0.300 

-0.274 

0.987 

0.991 

20 

B 

D/C 

W/w 

9.666 

8.578 

-0.301 

-0.268 

0.978 

0.963 

8 

C 

D/C 

W/w 

5.572 

4.263 

-0.389 

-0.328 

0.921 

0.965 

10 

D/C  =  drying/cooling  for  SWC/SFC  data 
W/w  =  wetting/warming  for  SWC/SFC  data 
=  (or/porosity)'1''® 


8  (*C) 


-0.0004  -0.004  -0.04  -0.4  -4.0  -40.0 


a.  Drying  curve  data  from  SWC  of  Iskandar  et  ai. 
(1979).  Cooling  curve  data  from  SFC,  qj  =  1.38  g/cm'. 
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b.  Wetting  curve  data  from  SWC  of  Iskandar  et  al. 
(1979).  Warming  curve  data  from  SFC,  =  1.38  g/cm'. 


Figure  l.  Semilog  plot  of  Wv  vs  <baw  and  6  for  Windsor  soil 
A  horizon. 
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sides  and  complete  inidal  saturation  should  be 
used.  The  two  cases  for  which  the  dry  bulk  densi¬ 
ties  (Cd)  are  most  similar  are  the  warming  and 
cooling  data  for  the  A  horizon  with  ed  =  1  -38  and 
the  warming  and  cooling  data  for  the  B  horizon 
with  cd  =  1.56.  In  both  cases,  the  high  degree  of 
saturarion  meets  the  air-free  requirement.  The  dry 
bulk  densities  used  for  the  unfrozen  water  content 
measurements  in  the  C  horizon  were  all  greater 
than  the  density  reported  by  Iskandar  et  al.  (1979) 
of  1.534.  The  specimen  that  came  closest,  ed  = 
1.79  and  85.8%  saturation,  was  used. 

Since  the  clay  content  of  the  Windsor  soil  is  less 
than  1%  (Iskandar  et  al.  1979),  the  transforma¬ 
tion  between  SWC  and  SFC  data  is  obtained  with 
eq  T2  of  Table  1.  Figures  1-3  are  plots  of  trans¬ 
formed  SFC  and  SWC  data  along  with  the  best  fit 


regression  to  eq  9.  Table  3  contains  a  summary 
analysis  for  each  plot.  SWC  data  that  appeared  to 
be  above  the  air-entry  value  (<£,w  <  <^)  were  ex¬ 
cluded  in  the  regression. 

Inspection  of  Figures  1-3  and  the  high  ri  values 
( >  0.9)  for  the  regression  fit  clearly  indicates  the 
success  of  eq  9  in  describing  both  SWC  and  SFC 
data  for  this  soil  and  the  suitability  of  eq  T1  of 
Table  1  for  relating  the  two  sets  of  soil  data.  The 
worst  case  (r*  =  0.921)  is  the  drying/cooling  data 
for  the  C  horizon  (Fig.  3a).  The  data  sets  were  not 
at  the  same  bulk  densities,  which  violates  the  re¬ 
quirement  for  comparison.  Additional  efforts 
failed  to  produce  a  soil  specimen  with  as  low  a 
density  (ed  =  1 .534)  as  reported  by  Iskandar  et  al. 
(1979). 

The  analysis  presented  in  this  report  extends  the 
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a.  Drying  curve  data  from  SWC  of  Iskandar  et  al. 
(1979).  Cooling  curve  data  from  SFC,  qj  =  1.56  g/cm'. 
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b.  Wetting  curve  data  from  SWC  of  Iskandar  et  al. 
(1979).  Warming  curve  data  from  SFC,  ea  —  1.56  g/cm'. 


Figure  2.  Semilog  plot  of  Wv  vs  <t>aw  and  ef°r  Windsor  soil 
B  horizon. 
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a.  Drying  curve  data  from  SWC  of  Iskandar  et  al. 
(1979).  Cooling  curve  data  from  SFC,  ea  =  1-79  g/cm ’. 
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b.  Wetting  curve  data  from  SWC  of  Iskandar  et  al. 
(1979).  Warming  curve  data  from  SFC,  =  1.79  g/cm'. 

Figure  3.  Semilog  plot  of  W„  vs  <baw  and  d  for  Windsor  soil 
C  horizon. 


range  of  the  SFC  data  to  higher  temperatures  than  CONCLUSIONS 
are  easily  obtained  with  present  techniques.  Corre¬ 
spondingly,  one  may  extend  the  range  of  the  SWC  In  this  report,  data  for  the  unfrozen  water  con- 

data  to  a  suction  of  greater  than  60  MPa.  In  addi-  tent  as  a  function  of  temperature  (SFC)  were 

tion,  the  new  interpretation  of  the  free  parameters  measured  and  analyzed  in  terms  of  a  power  curve 

in  the  power  curve  relationship  allows  determina-  relationship.  These  data  were  then  related  to  ice- 
tion  of  the  air-entry  value  from  air-free  SFC  data.  free  water  content  as  a  function  of  matric  suction 

It  is  of  interest  to  note  that  the  values  for  a  and  data  (SWC),  collected  nearly  a  decade  earlier, 

0  determined  from  the  SFC  data  alone  are  not  ap-  through  the  modified  Clapeyron  equation  and  sur- 

preciably  different  from  the  values  determined  face  tension  adjustments.  It  was  shown  that  when 

from  the  combined  SWC  and  SFC  data.  This  indi-  the  soils  were  strictly  similar  (i.e.  when  the  same 
cates  that  SFC  data  alone  are  good  indicators  of  bulk  density  was  used  for  the  SWC  and  SFC  data), 

the  behavior  of  this  type  of  granular  soil  and  may  then  a  single  relationship  between  moisture  con- 

be  employed  where  SWC  data  do  not  exist.  tent  and  the  ^variable  obtained  from  SWC  data 


and  air-free  SFC  data  described  these  data, 
demonstrating  the  validity  of  the  procedure  for 
granular,  colloid-free  soil. 

The  method  can  expand  the  utility  of  data  bases 
for  soil  physical  properties.  When  data  are  care¬ 
fully  collected,  and  experimental  conditions  duly 
noted,  the  results  compiled  for  one  intention  may 
be  used  by  others  for  different  objectives  at  later 
times. 
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APPENDIX  A:  SOIL  FREEZING  CURVE  DATA 


Table  Al.  Regression  parameters  to  eq  9  (i.e.  Wm  -  a(P)  for  unfrozen  water  content  data 
listed  below  on  a  percent  weight  basis  as  determined  by  NMR  and  initial  water  contents  and 
balk  densities. 


a 

a 

2 

r 

wT 

pd 

W_ 

vT 

e 

A  v 

c 

4.081 

4.758 

-0.253 

-0.307 

0.997 

0.993 

9.98 

1.63 

16.24 

38.0 

V 

c 

4.015 

4.608 

-0.259 

-0.300 

0.993 

0.993 

15.25 

1.62 

24.68 

38.4 

w 

c 

4.014 

4.814 

-0.277 

-0.326 

0.996 

0.979 

26.34 

1.48 

38.98 

43.7 

V 

c 

3.568 

4.155 

-0.282 

-0.317 

0.984 

0.986 

28.65 

1.54 

44.10 

41.4 

w 

c 

3.653 

4.130 

-0.281 

-0.309 

0.983 

0.996 

32.60 

1.38 

44.85 

47.5 

B  w 

c 

2.523 

2.697 

-0.322 

-0.325 

0.990 

0.984 

6.65 

1.74 

11.59 

35.3 

w 

c 

2.185 

2.357 

-0.291 

-0.311 

0.995 

0.970 

11.77 

1.78 

21.00 

33.8 

w 

c 

2.100 

2.375 

-0.304 

-0.331 

0.972 

0.997 

15.45 

1.79 

27.64 

33.5 

w 

c 

2.789 

2.911 

-0.354 

-0.348 

0.982 

0.980 

24.58 

1.55 

38.11 

42.4 

w 

c 

2.597 

2.564 

-0.364 

-0.356 

0.965 

0.980 

25.28 

1.56 

39.42 

42.0 

C  w 

c 

0.892 

0.983 

-0.357 

-0.405 

0.938 

0.987 

8.98 

1.75 

15.70 

35.9 

w 

c 

0.872 

0.897 

-0.393 

-0.391 

0.963 

0.973 

11.25 

1.80 

20.27 

34.1 

w 

c 

0.881 

0.853 

-0.396 

-0.399 

0.989 

0.988 

13.24 

1.79 

23.71 

34.4 

w 

c 

0.925 

0.989 

-0.384 

-0.419 

0.920 

0.984 

16.54 

1.79 

29.54 

33.4 

S 

42.7 
64.3 

89.1 

106.4 

94.4 

32.8 

62.1 
82.6 

89.9 

93.8 

44.7 

59.5 

68.9 

85.8 


* 


-  total  gravimetric  water  content 
(mass  water/mass  soil)*100 

-  total  volumetric  water  content 
(volume  water/total  volume) *100 

-  dry  density  (mass  soll/total  volume)  (g/cc) 

-  porosity  (volume  pores/total  volume)*100 

-  relative  degree  of  saturation  (V ^/e)*100 

-  warming  curve 

-  cooling  curve 
•  soil  horlson 

-  used  In  SVC/SFC  comparison 
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Data:  12/31/87  Tlaa:  16:03:57.69  Oparator:  PBB  Data:  12/31/87  Tlaa:  16:07:11.96 
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Data:  01/02/88  Tlaa:  14:49:35.82  Operator:  PBB  Date:  01/02/8B  Tiae:  14:50:40.63  Operator:  PBB 
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:  01/02/a*  Tlaa:  15:14:14.52  Operator:  PBB  Dace  01/02/8*  Tlaa:  15:18:42.89  Operator:  PBB 
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Data:  01/02/8*  Tlaa:  15:43:36.59  Oparator:  PBB  Date:  01/02/88  Tlae:  15:46:22.19  Operator:  PBB 
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D*t,:  01/04/88  XIm:  IS: 00: 38. 77  Operator:  PBB  Data:  01/04/88  Tina:  15:04:16.44  Oparator:  PBB 
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Data:  01/04/88  Tiae:  15:07:04.29  Operator:  PBB  Data:  01/04/88  Tlaa:  15:09:52.86  Operator:  PBB 
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Data:  01/04/88  Ti«a:  15:18:36.41  Operator:  PBB 
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APPENDIX  B:  ERROR  ANALYSIS 


Temperature 

Sample  temperatures  were  determined  from  a  thermocouple  immersed  in  the  constant 
temperature  bath  to  the  depth  of  the  soil  samples  and  measured  with  a  NESLAB  DR-2  digi¬ 
tal  readout  to  an  accuracy  of  ±0.1°C: 

A(Temperature)  =  ±0.1°C. 

Gravimetric  water  content 

Gravimetric  water  contents  were  determined  by  measuring  the  difference  between  wet 
sample  and  oven-dried  sample  mass: 

Wwet  ~  A^dry 
W  '  Mdry-AW 

where  A/wet  =  mass  of  wet  soil  and  tare 

A/dry  =  soil  mass  (oven-dried)  plus  tare 
Mare  *  mass  of  the  tare. 

The  error  involved  is  obtained  by  application  of  the  chain  rule: 
dtVw  dfVw 

=  aA^TAA/wet+  lM^AMdry+  'dM^eAM"Te' 

The  partial  differentials  are  found  to  be 


dW„ 

1 

dWwet 

Miry  —  Mare 

^W+l 

dMity 

Miry 

and 


dWv  _ 

W'w 

Mare 

Mdry  "Aftiit 

with 

AW, out  =  ±0.001  g 
AJtf|oj|  =  ±  0.001  g 
and 

AMajc  =  ±  0.001  g. 
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Unfrozen  water  content 

The  unfrozen  water  content  is  determined  from  the  ratio  of  measured  NMR  signal  to  the 
“theoretical”  NMR  signal  of  supercooled  water  at  that  temperature,  both  adjusted  for  back¬ 
ground  NMR  signal,  multiplied  by  the  gravimetric  water  water  content  of  the  ice-free  soil: 


UWC 


NMR-BKG 
TNMR-BKG  w 


where  UWC  is  the  calculated  unfrozen  water  content,  NMR  is  the  signal  measured  at  a  tem¬ 
perature,  TNMR  is  the  calculated  NMR  signal  for  a  theoretical  supercooled  soil  water  at  that 
temperature  obtained  from  a  linear  regression  of  the  above-0°C  NMR  data  and  BKG  is  the 
background  NMR  signal.  The  uncertainty  is  determined  by  applying  the  chain  rule: 


AUWC  = 


auwc 

3NMR 


ANMR  + 


3UWC 

3BKG 


ABKG  + 


auwc 

3TNMR 


ATNMR  + 


auwc 

aww 


AtVv. 


Solving  for  each  partial  gives 

auwc 

3NMR  =  TNMR-BKG 

aUWC  NMR -TNMR 
3BKG  *  (TNMR-BKG)1  w 

aUWC  _  NMR-BKG 
3TNMR  (TNMR-BKG)1  w 


aUWC  NMR-BKG 
dWm  ~  TNMR-BKG  ' 


The  uncertainty  of  each  term  associated  with  the  NMR  signal  is  found,  after  years  of  ob¬ 
servation,  to  be 


ANMR  =  ±2  counts 


ABKG  =  ±  2  counts 

and  because  of  the  near  perfect  linear  correlation  of  above  0°C  data,  the  uncertainty  of  the 
“theoretical”  supercooled  signal  is  assumed  to  be  the  same  as  the  others: 

ATNMR  ~  ±2  counts. 
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APPENDIX  C:  SOIL  WATER  CURVE  DATA 


Table  Cl.  Soil  water  curve  (SWC)  data  for  three  hori¬ 
zons  of  the  Windsor  soil  from  Iskander  et  ai.  (1979). 


Soil 

horizon 

Suction, 

*rm 

(bars) 

Volumetric 
water  content, 

w. 

Suction, 

4>«w 

(bars) 

Volumetric 
water  content, 

w. 

A 

0.000 

0.450 

0.588 

0.150 

0.025 

0.430 

0.294 

0.170 

0.049 

0.420 

0.132 

0.200 

0.103 

0.300 

0.049 

0.255 

0.216 

0.220 

0.000 

0.420 

0.348 

0.190 

0.S29 

0.170 

0.745 

0.150 

0.980 

0.140 

B 

0.000 

0.350 

0.980 

0.105 

0.01 5 

0.340 

0.588 

0.110 

0.025 

0.330 

0.402 

0.115 

0.049 

0.225 

0.1% 

0.135 

0.059 

0.205 

0.098 

0.150 

0.078 

0.185 

0.059 

0.165 

0.137 

0.135 

0.039 

0.180 

0.245 

0.150 

0.025 

0.200 

0.402 

0.135 

0.015 

0.240 

0.667 

0.120 

0.000 

0.335 

0.980 

0.105 

C 

0.000 

0.360 

0.980 

0.035 

0.020 

0.345 

0.662 

0.037 

0.054 

0.230 

0.392 

0.050 

0.078 

0.160 

0.1% 

0.070 

0.127 

0.125 

0.132 

0.090 

0.1% 

0.100 

0.064 

0.120 

0.392 

0.070 

0.000 

0.300 

0.529 

0.055 

0.667 

0.050 

0.980 

0.035 
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>*Vf.  /> 


A  facsimile  catalog  card  in  Library  of  Congress  MARC  format  is  repro¬ 
duced  below. 
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